Effects of optical and surface polar phonons on the optical conductivity of doped 

graphene 
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Using the Kubo linear response formalism, we study the effects of intrinsic graphene optical and 
surface polar phonons (SPPs) on the optical conductivity of doped graphene. We find that inelastic 
electron-phonon scattering contributes significantly to the phonon-assisted absorption in the optical 
gap. At room temperature, this midgap absorption can be as large as 20-25% of the universal ac 
conductivity for graphene on polar substrates (such as AI2O3 or Hf02) due to strong electron-SPP 
coupling. The midgap absorption, moreover, strongly depends on the substrates and doping levels 
used. With increasing temperature, the midgap absorption increases, while the Drude peak, on 
the other hand, becomes broader as inelastic electron-phonon scattering becomes more probable. 
Consequently, the Drude weight decreases with increasing temperature. 
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I. INTRODUCTION 

Since it was first isolated in 2004, 1 graphene, a material 
composed of a single layer of carbon atoms arranged in a 
two-dimensional (2D) honeycomb lattice, has attracted 
immense interest 2-4 due to its excellent transport and 
optical properties, 2,4-9 which make it an attractive can- 
didate for possible applications in nanoscale electronics 
and optoelectronics. 10-12 One particular field which has 
received considerable attention, both experimentally 13-16 
as well as theoretically, 17-25 is the optical (or ac) conduc- 
tivity in graphene, that is, the frequency-dependent con- 
ductivity. The main feature that can be observed in the 
optical conductivity is that for frequencies larger than 
twice the absolute value of the chemical potential fj,, the 



optical conductivity is roughly given by <7q 



7(4ft), 



the so-called universal ac conductivity. 13 ' 14 For frequen- 
cies below 2|//|, the optical conductivity is greatly re- 
duced, which can be explained within a single-particle 
model where transitions induced by photons with ener- 
gies hui < 2 1 I are forbidden due to Pauli's exclusion 
principle. The ability to tune the optical properties of 
graphene has been explored for use in broadband light 
modulators. 26-28 One figure of merit for this application 
is the modulation depth of the optical absorption. In 
experiments, however, one does not observe the optical 
conductivity to vanish completely, as one would expect 
from the simple single-particle argument given above. In 
addition, a substantial Drude weight loss has been re- 
ported in graphene on Si02- 16 

To describe this behavior, mechanisms involving dis- 
order and/or phonons, both of which can account for 
a finite absorption below 2|/x|, have been studied the- 
oretically in both monolayer 20-22,25 ' 29 and bilayer 30 ' 31 
graphene. In addition to these single-electron effects, 
excitonic effects 23 as well as effects arising from the 
Coulomb interaction 24 have also been considered, but 
were found to have a negligible effect on the midgap ab- 



sorption in heavily doped samples. Moreover, the op- 
tical conductivity in the presence of a magnetic field, 
the so-called magneto-optical conductivity, has also been 
investigated theoretically, 32,33 with Ref. 33 taking into 
account the coupling between electrons and Einstein 
phonons. 



Besides the aforementioned studies on the optical 
conductivity, the role played by different phonons has 
also been studied in the context of heat dissipa- 
tion mechanisms 34-36 and current/velocity saturation in 
graphene, 37 which plays an important role in electronic 
RF applications 11 and also for transport 38,39 in the sim- 
ilar system of carbon nanotubes. Inelastic scattering ei- 
ther by intrinsic graphene optical phonons 40 or surface 
polar phonons 37,41-45 (SPPs) is thought to give rise to 
the saturation of the current in graphene and affect the 
low field carrier mobility. 46,47 However, from transport 
experiments alone it is difficult to identify the role played 
by SPPs from the polar substrates because of the compli- 
cations arising from charge traps which can be populated 
thermally 48 or by the high electrical fields. 49 



Here, we show that the temperature dependence of the 
midgap absorption can be significantly stronger in the 
presence of SPPs as compared to suspended graphene or 
graphene on a non-polar substrate such as diamond-like 
carbon. Our main goal in this manuscript is to study the 
optical conductivity in the presence of phonons. While 
the impact of optical phonons has been studied in sev- 
eral earlier works, 20-22 the effect of SPPs on the optical 
conductivity in graphene has yet to be analyzed. In this 
paper, we use linear response theory to derive a Kubo for- 
mula for the optical conductivity, which is then evaluated 
for suspended graphene as well as graphene on different 
polar substrates, where SPPs arc present. 
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II. MODEL 

To describe the electronic (single-particle) band struc- 
ture of graphcne, we use the Dirac-cone approximation, 
where the Hamiltonian can be written as 



Ho 



E 



(1) 



with ek = hv F k. Here, k, s, and v denote the momen- 
tum, spin, and valley quantum numbers, A the conduc- 
tion (A = +1) and valence (A = —1) bands, c\ ksv and 
cxksu the corresponding creation and annihilation opera- 
tors, and v F « 10 8 cm/s the Fermi velocity in graphene. 4 
Since the goal of this work is to study and compare 
the effects of several different phonons on the optical 
conductivity of graphene, we need to take into account 
the interaction with these phonons. A general phononic 
Hamiltonian reads as 



H 



= ^hw K (q)p qA p qA , 

q,A 



(2) 



where different phonon branches are labeled as A, the 
phonon momentum as q, and the corresponding frequen- 
cies and creation (annihilation) operators as co A (q) and 
PqA (PqA)- Whereas Eqs. (1) and (2) describe isolated 
systems of electrons and phonons, respectively, the cou- 
pling between those systems is given by 



tfc-ph = E E M™'',A(k,q) (£ qA +p q A 



Aksv A'qv'A 



(3) 



X ^A(k+q)st)^ A ' k ™'' 



where A (k, q) is the electron-phonon coupling ma- 
trix element. 50 Hence, the total Hamiltonian of our model 
reads as 



H = H c +H, 



pii 



-ph- 



(4) 



In this work, we investigate two different types of 
phonons that couple to the electrons in graphene: intrin- 
sic graphene optical phonons and SPPs, that is, surface 
phonons of polar substrates which interact with the elec- 
trons in graphene via the electric fields those phonons 
generate. Intrinsic graphene acoustic phonons are not 
included in our model because their effect on the optical 
conductivity is negligible as has been shown in Ref. 20. 
The dominant electron-optical phonon coupling 51-54 is 
due to longitudinal-optical (LO) and transverse-optical 
(TO) phonons at the T-point and TO phonon at the 
i^-point. In the vicinity of the T point, the disper- 
sion of both, LO and TO phonons (denoted by r LO 
and r TO ), can be approximated by the constant energy 
hoj r rj 197 meV. Near the K and K' points, on the 
other hand, only the TO phonon (denoted by K TO ) con- 
tributes to the electron self-energy and its dispersion in 
this region can again be assumed as flat, huj K w 157 



meV. 



67 



Furthermore, 

AA 



we need to know the products 



M A = £M$ A (k - <l'<l) M iv<,A( k > -<J) of the electron- 

V 

phonon coupling matrix elements for each phonon branch 
in order to calculate the electron self-energy. The cou- 
pling of both r phonons is given by 51,52,54 



Mr 



2NM r .oj, 



(1 + AA')<U<, (5) 



where N is the number of unit cells, M c the carbon mass, 
and D T sa 11.2 eV/ A the strength of the electron-phonon 
coupling. 42 The coupling of the K TO phonon mode to the 
electrons in graphene is twice as large as that of phonons 
at the T-point 51-53 and is described by 68 



hDl 



2NM r u h 



1 + AA' - A (Ae- ie + A'e ie ) 



(6) 

where 6 = 9 k — # k _ q and 9 k = arctan (k x /k y ). Our choice 
of the optical phonon deformation potential lies between 
the LDA results of Rcfs. 51 and 52 and the GW results 
from Ref. 55. A smaller electron-intrinsic optical phonon 
coupling as suggested in Ref. 56 would further reduce the 
absorption below 2|^i| in suspended graphene. 

Here, we include SPPs in our model as follows: There 
are two surface optical (SO) phonons in polar substrates 
that interact with the electrons in graphene and whose 
dispersion can again be approximated by substrate- 
specific, constant frequencies lo SOi and w S02 , and their 
electron-phonon coupling matrix elements read as 



43.47 



M A = 



rVR 2 



NAq 



2qzn 



1 + AA' + A (Ac -19 + A'e ie ) 



5w, 
(7) 

where e = |e| is the absolute value of the electron charge, 
A = 3\/3a 2 /2 the area of the graphene unit cell, a « 1.42 
A the distance between two carbon atoms, zq w 3.5 A the 
van der Waals distance between the graphene sheet and 
the substrate, and the Frohlich coupling Fl(q) describes 
the magnitude of the polarization field. 57 The Frohlich 
coupling is given by 34,58 



and 



FioM 



2tt 



1 



1 



2rr 



£i + e(g) e + e(q) 



1 



(8) 



(9) 



with the optical, intermediate, and static permittivi- 
ties e x , e t , and e of the substrate as well as the 
static, low temperature dielectric function e(q) = 1 + 
27re 2 n g (g,w = 0) /(nq), where k is the background di- 
electric constant and n g (q, w) the polarization function 
of graphene as calculated in Refs. 59,60. The dielectric 
function e(q) accounts for the screening of the Coulomb 
interaction in the graphene sheet above the polar sub- 
strate. If the effect of screening is to be disregarded, we 
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Al 2 3 a 


h-BN 6 


Hf0 2 c 


SiC d 


Si0 2 e 


Sa 


12.53 


5.09 


22.0 


9.7 


3.90 


£i 


7.27 


4.575 


6.58 


- 


3.36 




3.20 


4.10 


5.03 


6.5 


2.40 


hcjso-t [meV] 


56.1 


101.7 


21.6 


116.0 


58.9 


fcso 2 [meV] 


110.1 


195.7 


54.2 




156.4 


*so n KV] 


0.420 


0.258 


0.304 


0.735 


0.237 


Flo, [meV] 


2.053 


0.520 


0.293 




1.612 



TABLE I: Optical, intermediate, and static permittivities as 
well as frequencies and bare Frohlich couplings for the SPP 
scattering on the substrates AI2O3, hexagonal BN, HTO2 SiC, 
and SiC>2. 
a Refs. 58,61 
b Refs. 42,61,62 
c Refs. 42,58,61 
d Refs. 42,63 

e Ref. 42, which uses averages of values from Refs. 43,47,58. 



set e(q) = 1 in Eqs. (8) and (9), for which we obtain the 
bare Frohlich couplings presented in Table I. 

Employing standard diagrammatic perturbation 
theory 50,64,65 and inserting the specific expressions for 
the matrix elements M^ A (k, q) (for more details, we 
refer to the Appendix A), we find that, up to first 
non-vanishing order, the electronic spectral function is 
diagonal in the four quantum numbers A, k, s, and v 
and is given by 



Ax (k,w) 



- 21m 



1 



u + i0+ - [Ae k - fi + £ A 0, lj + i0+)] /K 



(10) 

where £ A (k,iv n ) denotes the imaginary-time self-energy 
at the imaginary (fcrmionic) frequency \v n . In the lowest 
order, the contribution to the self-energy due to phonons 
is just the sum of the contributions from the different 
phonons A of our model. Moreover, we include scattering 
at Coulomb impurities in our model by adding the con- 
tribution S Co (k) = —ih/[2T(k)] to the self-energy, where 
we use the transport scattering time r(k) as calculated 
in Ref. 66, and the total self-energy reads as 



£ A (k, w n ) -E c ° (*) + S r (nv„) + £ K (w n ) 



(k,iv n ) 



(k,iu n ) 



(11) 



have the form 
Im [S A (u> + 



i0 + )] = 

3 (hw A — hu>)] g A (hui +11 — hui A ) 
i FD (hu) A + huj)] g A (tuv + /J + hw A ) 



(12) 



with the functions g A (e) = — AD 2 \e\/(2M c hu} A vl). 



Similarly, the effect of the 
SOi, SO2) is described by 



two SPP modes (A 



Irnpt (k, 



where 



K(k,e) 



fi0 + )] = 

(huj A — fno)] (k, fvjj + /x — hu A ) 
; FD (hio A + hu))] h'l (k, huj + + huj A ) . 



(13) 



2 71 



7re 



2(hv F ) 2 



dO 



Fl(q)e- 2qzo 



F Ae cos 9) 
(14) 

with q = yj( 2 + t't- 2ee k cos 8/(hv F ). In Eqs. (12) 
and (13), we have introduced the Fermi-Dirac and Bosc- 
Einstein distribution functions, n FD/BE (e) = l/[exp( J 8e)± 
1], where /3 = l/(k B T) (with T and fc B being the temper- 
ature and the Boltzmann constant, respectively), n A = 
n BE (fno A ), and the chemical potential /x = fi(T). In 
the following, we ignore the effect of polaronic shifts, 
that is, the real parts of the self-energies and set 
S A (fc,u; + i0 + ) ee iIm[S A (fc,a; + i0 + )]. The total (re- 
tarded) self-energy in our model is thus completely imag- 
inary. 

As detailed in the Appendix B, the real part of the 
optical conductivity can be calculated from the spectral 
function (10) via the formula 69 



a v 



E 



da/ 



dkkAx (k,ui') 
x Ay (k, oj' + oj) [n FD (fvJJ 1 ) — n FD (Huj ■ 



fa')], 

(15) 

which includes the universal ac conductivity ao = 
e 2 /(4fi). In the following, we will use Eqs. (10)-(14) to 
numerically calculate the spectral function, which is in 
turn used to calculate the real part of the optical con- 
ductivity numerically via Eq. (15). Those calculations 
are conducted for several different substrates (as well as 
suspended graphene) with the corresponding parameters 
summarized in Table I. 



The contribution S Co (k) has been added to the self- 
energy to model the lineshape of the Drude absorption 
peak. Throughout this work we use the impurity con- 
centration of TT-j = 5 x 10 11 cm -2 , which is low enough 
to not affect the midgap absorption of graphene on polar 
substrates significantly. 

The imaginary parts of the contributions from the op- 
tical phonons at the T and K points (A = T, K) to the 
retarded self-energy depend only on the frequency and 



III. RESULTS 

In Fig. 1, the optical conductivities 70 (for a fixed 
chemical potential fx = 0.3 eV) at two temperatures 
T = 1 K (inset) and T = 300 K are shown for sus- 
pended graphene as well as graphene on several differ- 
ent substrates: AI2O3, hexagonal BN, HfC>2, SiC, and 
Si02- For comparison, we have also included the optical 
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jl 1 \ Drude peak 
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FIG. 1: (Color online) Calculated frequency dependence 
of the (real part of the) optical conductivity of suspended 
graphene and graphene on several different substrates for 
rii = 5 x 10 11 cm~ 2 and fj, = 0.3 eV at T = 300 K and 
T= 1 K (inset). 
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FIG. 2: (Color online) Calculated frequency dependence of 
the (real part of the) optical conductivity of graphene on a 
Si02 substrate for rii = 5 x 10 11 cm~ 2 and several different 
chemical potentials at (a) T = 1 K and (b) T = 300 K. 



conductivity of suspended graphene and graphene with- 
out any phonon contribution (at k = 1). Figure 1 has 
been calculated using the parameters given from Table I, 
the dielectric function e(q), and k = (1 + e )/2 as the 
background 60,66 dielectric constant. 

The profiles in Fig. 1 illustrate the main features that 
the effect electron-phonon coupling has on the optical 
conductivity: Whereas there is a gap with a width 2\/j.\ 
in the absorption spectrum of the purely electronic single- 
particle model, where direct transitions between the elec- 
tronic states in the conduction and valence bands are for- 
bidden for energies < ftw < 2\/i\ due to Pauli blocking, 
there is a finite absorption in this region in the presence of 
phonons. This finite absorption is largely due to pho non- 
assisted transitions which give rise to distinct sidebands, 
the onsets of which can clearly be distinguished from the 
Drude peak at low temperatures and low impurity densi- 
ties (see the inset in Fig. 1). If the photon energy exceeds 
2\fi\, direct (interband) transitions become possible, re- 
sulting in a steep rise of the optical conductivity. 

At higher temperatures, one can see that the Drude 
peak is broadened as more phonons become available and 
electron-phonon scattering becomes more probable. For 
Hf02 and AI2O3 substrates, the phonon sidebands merge 
with the Drude peak resulting in a very broad Drude 
peak at room temperature. Furthermore, the profiles of 
the optical conductivity are much smoother compared to 
those at T = 1 K and distinct onsets of phonon side- 
bands can no longer be observed as the profiles of the 
optical conductivity are smeared out by thermal broad- 
ening. Finally, Fig. 1 shows that the so-called "midgap 
absorption" , that is, the absorption at fko = fi, is sig- 
nificantly enhanced for graphene on polar substrates as 
compared to suspended graphene or graphene on nonpo- 
lar substrates: Whereas the midgap absorption at room 
temperature is about 5-6% of Co for suspended graphene, 
it can be as high as 20-25% of ao for graphene on HTO2 
or AI2O3. Hence, the midgap absorption strongly de- 



pends on the particular polar substrate used and is, in 
particular, determined by the interplay between the SPP 
frequencies w SOi and the Frohlich couplings F<? : the 
smaller the w SOi or the larger F s 2 ., the larger is the 
midgap absorption. 

Impurity scattering has also an influence on the 
midgap absorption: By calculating the absorption spec- 
tra using just the Coulomb impurity scattering with 
K = 1 and no phonons for /x = 0.3 eV the Kubo model 
shows that while at T = 300 K the midgap absorption 
is 2.3% for rii — 5 x 10 11 enr 2 (see also Fig. 1), it in- 
creases to 4.3% for n, = 10 12 cm" 2 . At T = 500 K, 
these numbers are 5% and 7%, respectively. We note 
that an obvious first estimate of the midgap absorp- 
tion could have been obtained by using a Drude model 
a(u))/ao = 4|/i|/[7rfir(o; 2 + 1/t 2 )], where r is a scat- 
tering time. For n-, = 5 x 10 11 cm~ 2 and ?i ; = 10 12 
cm" 2 , the Coulomb scattering mobilities arc 7 w 13000 
cm 2 /(Vs) and 7 ~ 6500 cm 2 /(Vs), respectively, corre- 
sponding to midgap absorptions from the Drude model 
of <j(^/K)/cr G fn 4hevp/(7r^fi 2 ) = 0.7% and 1.4%, respec- 
tively. Thus, this simple estimate using the Drude model 
significantly underestimates the results obtained from the 
full calculations. Indeed, the deviations between the esti- 
mate from the Drude model and the full Kubo formalism 
calculation become even more pronounced if phonons are 
taken into account. 

Figure 2 shows the optical conductivity for graphene 
on a SiC>2 substrate at different temperatures and chem- 
ical potentials. Apart from the trends in the behavior of 
the optical conductivity discussed above, one can clearly 
see different gaps in the absorption spectrum, given by 
2\/j,\ for each chemical potential. Another feature that 
can be discerned from Fig. 2 is that the maximal value 
of the phonon-mediated absorption in the gap increases 
with increasing chemical potential (doping level). More- 
over, we note that due to the electron-hole symmetry of 
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FIG. 3: (Color online) Calculated dependence of the midgap 
absorption (at hwo = /j) on the chemical potential for sus- 
pended graphene and graphene on several different substrates, 
m = 5 x 10 11 cm -2 , and T = 300 K. 



FIG. 5: (Color online) Calculated temperature dependence of 
the midgap absorption (at hwo = fi) for suspended graphene 
and graphene on several different substrates, rii — 5 x 10 11 
cm -2 , and /j = 0.3 eV. 
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FIG. 4: (Color online) Calculated frequency dependence 
of the (real part of the) optical conductivity of suspended 
graphene and graphene on several different substrates for 
Hi = 5 x fO 11 cm" 2 and fj, = 0.3 eV at T = 300 K and 
T = 1 K (inset) if bare Frohlich couplings are used. 



the Dirac Hamiltonian, the profiles of the optical conduc- 
tivity would look the same for p-doped graphene. The 
dependence of the midgap absorption on the chemical 
potential at room temperature is shown in Fig. 3, again 
for suspended graphene and graphene on several different 
substrates. In the region studied here between /i = 0.2 
eV and /i = 0.4 eV, the midgap absorption decreases with 
increasing chemical potential for graphene on substrates, 
with the decay being most pronounced for Hf02- The dif- 
ferences in the shape of the phonon mediated gap absorp- 
tion reflect the momentum dependence of the clcctron- 
phonon interaction in Eq. (5)-(7). 

In particular, Fig. 4 reveals a striking difference in the 
absorption if we use a bare unscreened Frohlich couplings. 
If screening is not accounted for, we find that the opti- 
cal conductivity in the optical gap is greatly enhanced 
compared to the situation where screening is used. The 
most noticeable feature if the bare Frohlich coupling is 



used, is that a second clearly distinguishable phonon side- 
band peak (due to the SPPs) can now be observed in the 
absorption spectra even at room temperature for SiC>2 
and SiC substrates. For BN substrates, one can even 
find two such peaks at room temperature. We suggest, 
therefore, that measurements of the midgap absorption 
in graphene on different substrates could help to clarify 
the still controversial 34,43 ' 44 ' 58 ' 61 issue concerning the ef- 
fect of screening on the electron-SPP coupling strength. 

We also investigate the temperature dependence of the 
midgap conductivity (that is, at Hujq = fi) within our 
model in Fig. 5 on several different substrates. At tem- 
peratures below 100 K, the midgap absorption does not 
depend strongly on the temperature. At about 100 K, 
an increase of the optical conductivity at hwo = /i is pre- 
dicted to take place. Also, the smaller the energy of the 
dominant phonon contributing to the gap absorption, the 
stronger is the temperature dependence in Fig. 5. 

Finally, we relate the midgap absorption to the spec- 
tral weight of the Drude peak. Describing the graphene 
optical conductivity in the non-interacting single-particle 
picture, the spectral weight of the bare Drude peak is 

Iq = J duja(uj) — 2\fi\<Jo/h = Do/2, where Do is the 
o 

bare Drude weight and uj' is some characteristic fre- 
quency much larger than the scattering rate, but smaller 
than both the lowest energy of the optical phonon fiw opt 
and 2\fi\. In the presence of phonons, the total spec- 
tral weight has to be conserved. The spectral weight 
contribution due to the midgap absorption can be ap- 

proximated at low temperatures as J gap = J dwa(uj) ~ 

"opt 

a<7o(2|/i| — huj opt )/h, where a is the averaged value of 
cr((jj), that is, of the real part of the optical conductiv- 
ity, in units of uq- If wc further assume that the entire 
spectral weight lost at the Drude peak is transferred to 
the optical gap and that Huj opt <C 2\[i\, within this picture 



6 





20- 




15- 


o 

c 


10- 


3= 
C 


5- 




ALO, substrate 

. _ . . . SiO,, substrate 

suspended — " 




Q 80- 

a 

70-^ 

20 ; ( C ) 

> 15- 
fl 10- 
S 5 _ 

o 



0.2 



0.4 



0.6 



0.8 



FIG. 6: (Color online) Calculated dependence of the (a) ab- 
sorption at htoo = /J., (b) the Drude weight D, and (c) the 
inverse scattering time h/r on the relative strength of the 
electron-phonon coupling 7 rc i for suspended graphene and 
graphene on A1 2 3 and Si0 2 , n t = 5 x 10 11 cm' 2 , T = 300 
K, and \i = 0.3 eV. 



the remaining Drude weight D can then be estimated as 
D I Do = {Iq — I g a,p)/Io ~ 1 — ol. Thus, from this con- 
sideration we expect that, as a increases with increasing 
temperature or strength of the electron-phonon coupling, 
the Drude weight is reduced. 

Figure 6 shows (a) the absorption in the optical gap 
as well as the fitted (b) Drude weight and (c) inverse 
scattering time for suspended graphene and graphene on 
A1 2 3 and Si0 2 substrates with n = 0.3 eV, T = 300 
K, and 7ii = 5 X 10 11 cm~ 2 as functions of the relative 
strength 7 reI of the electron-phonon coupling. Here, the 
optical conductivities have been calculated by scaling the 
(products of) electron-phonon coupling matrix elements 
with Ma — > 7 rcl AfA, and the Drude weight D as well 
as the inverse scattering time 1/r have been extracted 
from the optical conductivity by fitting the Drude peak 
to a Lorentzian. As expected from the argument given 
above, the Drude weight decreases with increasing 7 rel , 
although the simple relationship between reduced Drude 
weight loss and the midgap absorption as 1 — a does not 
hold. This is because the midgap absorption does not 
coincide with the averaged value of <j(oj) in the gap and 
because some of the spectral weight from the Drude peak 
is transferred not only to the optical gap, but also to the 
spectral region Hu> > 2\fi\. 

With increasing 7, el , electron-phonon scattering be- 
comes more probable and consequently the scattering 
time decreases as can be seen in Fig. 6 (c). The corre- 
sponding increase of 1 /r is most pronounced for graphene 
on AI2O3 and least pronounced for suspended graphene. 
Finally, we remind the reader that for ?i; = 5 x I0 11 enr 2 
impurity scattering also contributes to the absorption in 
the optical gap, which can be discerned from the finite 



absorption at 7 rel = 0. Because suspended graphene and 
graphene on different substrates each possess different 
background dielectric constants and thus different trans- 
port scattering times, the residual values at 7 rcl = are 
also different. 



IV. CONCLUSIONS 

We have studied the effects of intrinsic graphene op- 
tical phonons and SPPs on the optical conductivity of 
doped graphene. Our focus has been on the absorption at 
frequencies Huj < 2\[i\, where optical transitions are for- 
bidden due to Pauli blocking in a clean system (at T = 0), 
but which can occur if phonons are present, giving rise 
to phonon sidebands. Here, we have found that inelas- 
tic phonon scattering contributes significantly to the ab- 
sorption in the optical gap and strongly depends on the 
substrate used: At room temperature (and /j, = 0.3 eV), 
the midgap absorption, which is mainly due to intrinsic 
optical phonons, amounts to about 5-6% of the universal 
ac conductivity for suspended graphene or graphene on 
non-polar substrates, while the midgap absorption can 
be as large as 20-25% of <jq for graphene on polar sub- 
strates (such as AI2O3 or HfC^) due to the smaller SPP 
energy and strong electron-SPP coupling. Moreover, the 
midgap absorption depends on the doping level and de- 
creases with increasing \fj,\, while the maximal value of 
the sideband absorption at low temperatures increases. 
We have also investigated the temperature dependence 
of the midgap absorption which increases with increas- 
ing temperature. The Drude peak, on the other hand, 
becomes broader with increasing temperature as inelas- 
tic electron-phonon scattering becomes more important. 
Consequently, the Drude weight decreases with increas- 
ing temperature due to the stronger phonon coupling. 
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Appendix A: Self-energy and Green's Function 

We use standard diagrammatic perturbation theory 
(with the unperturbed Hamiltonian H c + H p h and the 

perturbation 77 c _ p h) to calculate the electronic Matsub- 
ara Green's function of the system described by Eq. (4), 

hp 



Q% (k>„) = - / dr(T CAk™(T)cI, k8 „(0) 



(Al) 
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where r and \u n denote the imaginary time and 
(fermionic) frequency, (...) the thermal average, T 
the imaginary time-ordering operator, and ft — 
l/(fc B T). 50 ' 64,65 By solving the corresponding Dyson 
equation, we can express the electronic Green's func- 
tion via the self-energy Ti^'.w' (k,if n )- In writing down 
Eq. (Al), we have used that, due to the conservation of 
momentum and spin, the self-energy and thus the Green's 
function are diagonal in k and s and do not depend on s 
due to spin degeneracy. 

Up to the first non-vanishing order (and omitting the 
tadpole diagram, which yields a purely real self-energy 
contribution that can be absorbed in the chemical po- 
tential) , we find the electronic self-energy due to phonons 



for arbitrary matrix elements M^, A (k, q) to b 



.,50 



(A2) 



„, (k>,0«£i 2 A L/ (k>n) 

= \ XJ ^A(k-q,q)il4 A / A (k,-q) 

A,q,A,{) 

HqA + 1 - /Fp(Ae|k-q|) 

w n - (Ae| k _ q | - n)/H- w A (q) 

» qA + /F D (A6| k -q|) 

iv n - (Ae| k _ q | - /j,)/h + u} A (q) 



with the Fermi-Dirac and Bosc-Einstcin distribution 



functions, / FD (e) 



/z) and n q 



[haj A (q)}, 



and the chemical potential /i = /i(T) at the temperature 
T. Thus, in the lowest order, the self-energy is simply 
the sum of the contributions from the different phonons 
A. 

In order to calculate the self-energy, we need to 
know the products of the electron-phonon coupling 
matrix elements entering Eq. (A2), S^VcaC^ — 

V 

q,q)Mj$' A (k, -q), for the T, K, SOi, and S0 2 modes. 
Since ra q r TO = ?^ q r LO , it follows from Eqs. (A2) and (5) 
that the contribution due to the optical phonons near 
the r point is diagonal in the valley and band quantum 
numbers. We calculate the self-energy by using the trans- 
formation k' = k — q and replacing the sum ^ by the 

k' 

2D integral S/(2tt) 2 J dk'k' J d(9, where S = NA and 9 is 
chosen to be the angle between k and k'. 

Equation (6) describes the coupling of the K TO mode 
to the electrons in graphene. As above, the self-energy 
contribution from the K phonons is calculated by in- 
troducing k' and writing the sum as a 2D integral. 
After performing the integration over the angle 9, the 
terms containing e ±l9 in Eq. (6) vanish and consequently 
the contribution to the self-energy is diagonal with re- 
spect to the band quantum number. Moreover, Eqs. (6) 
and (A2) make it clear that, even though the matrix ele- 
ment M™i Kt (k, q) describes intervalley scattering, the 
second-order contribution from the K phonon to the self- 
energy is diagonal also in the valley quantum number. 



Using the SPP coupling matrix elements (7), the con- 
tribution from each SPP to the self-energy is again cal- 
culated by introducing k' and writing the sum as a 2D 
integral. Then, the angular integration for the offdiago- 
nal elements with respect to the band indices A and A' 

is of the type J d9 sin9f(cos9), where the function / 

o 

depends only on cos# and on the SPP considered, and 
vanishes for screened as well as unscreened SPPs. Thus, 
the lowest order contribution from each SPP to the self- 
energy is also diagonal in the band and valley quantum 
numbers. 

Combining the results discussed so far, the total con- 
tribution from all phonons is given by 



Ef , (k, w n ) = Waa'E? (k, w n ) , 

where Ef (k, \v n ) is given by 

Ef (fc,hv„) =E r K) + E K K) 

+ Ef! (fc>„) + £f 2 (fc> n ), 



(A3) 



(A4) 



and each individual contribution is calculated from 
Eq. (A2) as described above for A = A' and v = v'. Here, 
we find that, in contrast to the contributions E^ Ql (fc, iv n ) 
and E^° 2 (k,iv n ), the contributions from the graphene 
optical phonons, E r (iv n ) and E K (ii/„), do not depend 
on the band or momentum k. 

The contribution due to Coulomb impurity scattering 
reads as 



where 
E c °(fc) 



-m 



2r(fc) 



17T7li 



d 2 fc' 

(2tt)= 



2ne 2 



nqe(q) 



(A5) 



S(ek - £fc') 



(1 - cos0)(l + cos( 



(A6) 

with 9 = 9 k — 8 k > , q = k — k', the dielectric function e(q) , 
and the impurity concentration tv 66 Since this contribu- 
tion is also diagonal, the total self-energy 



(k,kv„) = 5 VV >5. 



ih 



2r(fc) 
= (5^'^aa'Ea (k,iu n ) 

and, consequently, the Green's function 



(A7) 



iv>„ - [Ae k - m + E A (fc, u> n )] /h (A8) 



are diagonal with respect to A and v in our model. Fi- 
nally, the spectral function is obtained from the Green's 
function via 



Ax (k,u) = -2lm[G x (fc,w + i0+)] 



(A9) 
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In this work, we are interested only in the imaginary 
parts of the retarded self-energy. Upon replacing iu n by 
w + iO + in Eq. (A2), the imaginary part of each contribu- 
tion A in Eq. (A2) contains a Dirac-<5 function [since there 
is no contribution from Im(MA.) as discussed above]. Af- 
ter introducing k' and writing the sum as a 2D integral, 
the Dirac-5 function can be used to calculate the k 1 inte- 
gral, which then yields Eq. (12) for the T and K phonons 
and Eqs. (13) and (14) for the SOi and SO2 phonons. 



Appendix B: Kubo formula for the optical 
conductivity 

1. Current density operator 

Our starting point in the derivation of a Kubo for- 
mula for the optical conductivity is the current operator. 
In the presence of an arbitrary magnetic vector poten- 
tial A(r), the (first-quantized) 2D Dirac Hamiltonian of 
graphene reads as 4 



H c = U F 7.7T 



(Bl) 



with 7r = p + eA(r) being the 2D kinetic momentum op- 
erator, p the 2D momentum operator, and the matrices 
lx = &x <S> 1, ly = &y <S> r z , j z = 0, where 1 is the 2x2 
unity matrix and <x and r are Pauli matrices referring to 
the A/B sublattices and the K/K' points, respectively. 
As discussed in Ref. 4, the 2D momentum k and the 
valley K/K' are good quantum numbers and the Hamil- 
tonian (Bl) has the (valley-degenerate) eigenvalues 



£A(k) = Aek = Xhv F k 
and the corresponding eigenstates 

near the K point and 





(B2) 



(B3) 



(B4) 



near the K' point, where S ~ NA denotes the sur- 
face area of the graphene sample, A = ±1, 6 k = 
arctan (k x /k y ), and 



Xk 



(B5) 



For an arbitrary (normalized) state ^(r), the energy 
expectation value as a functional of the vector potential 
A(r) is given by 

£[A]=W d a rK(x,v)(&*) *fi(z,y), (B6) 



where the sums over a and (3 refer to the matrix 7. The 
charge current density j(r) of this state ^(r) can be de- 
termined by a variational method: 

SE = E [A + 6 A] - E [A] = - J d 2 r j(r)c5A(r), (B7) 

which yields 

j(r) = -ev F Y, KM h) a p ( B 8) 

a/3 

Promoting the wave functions in Eq. (B8) to field oper- 
ators, using the cigcnbasis given by Eqs. (B2)-(B4), and 
taking into account the spin degeneracy, the charge cur- 
rent density operator can be determined as 



j(q)= Y d lA'(k,q)cIk s ^A'(k+q), 



kAA'st? 

in reciprocal space and as 



(B9) 



(BIO) 



in real space. Here, the dipole matrix elements read as 



<y,»(k,q) = 
Cv, B (k,q) 



-ev w 
2 



y e i(e k +e k+q )/2 + ^ c -i(e k +e k+q )/2 

y e i(0 k +e k+q )/2 _ Ae -i(e k +e k+q )/2 

d^,. / ,(k J q)= [d* w ,(k,q)]*. 

(Bll) 



2. Kubo formula 



If an external electric field 



E(r,i) = - 



dA(r,t) 
dt 



(B12) 



is applied to the system considered here, its effect can be 
described by 



fiext(*) = - / d 2 rj(r).A(r,i) 



(B13) 



with the charge current density operator (BIO). The to- 
tal Hamiltonian of the problem then reads as H + H ext (t) , 
where H is given by Eq. (4). 

Using linear response theory [for the unperturbed 
Hamiltonian H and the perturbation H cxt (t)] and con- 
ducting a Fourier transformation with respect to the time 
and position, 50 ' 64,65 we find that the current density due 
to the external field is given by 

*<ja(q,w)> =-\YKpMMw) ( B14 ) 

p 
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with Il^g (q, u>) being the (Fourier transformed) retarded 
current-current correlation function and a and (3 refer- 
ring to the in-plane coordinates x and y. The retarded 
correlation function H^p (q, w) can be related to the 
imaginary-time correlation function 



where 



Z^'(k,q) = d*,, a (k,q) [< v ,,(k,q)]* 



C;,„(k,q) (k,q)] : 



(B18) 



U af3 (q,iw„) = -- 



dr(T 



Ja(q, r)5>(-q,0) 



by Kb (q. 



H a/ 3 (q, w 



(B15) 

i0 + ), that is, by replacing 



iw„ with w + iO+ in Eq. (B15). 50 - 64 ' 65 Here, iw„ denotes a 
bosonic frequency. Hence, the Kubo formula for the real 
part of the conductivity reads as 



is a real number. 

Here, we are interested in the response to a uniform 
field, that is, in the case q = 0, for which Eq. (B17) 
becomes 



Im 



Re [a a p (q, ui) 



nf,(q, 



(B16) 



If vertex corrections due to phonons in Eq. (B15) are 
ignored, the phonon-dressed Green's functions given by 
Eq. (A8) expressed via their spectral functions (A9), and 
the sum over k rewritten as a 2D integral, we arrive at 



Re [a aB (0, to)} = S a pa (u) 



(B19) 



Re [cr Q(3 (q,w)] 



XX' 



d 2 fc 

(2tt) 2 



~2~7 



x £)**'(k,q) [n FD (fhu') - n FD (frw + ho')} 
x Ax (k,u')Ax> (\k + q\^T~ 



(B17) 



with a (u>) given by Eq. (15). In order to obtain Eqs. (15) 
and (B19), wc have used that only (k, 0) depends on 
the angle of the fc-intcgration in Eq. (B17) for q = and 

that / d0 k D**'(k,O) = ire 2 vj5 a p. 
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